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We have studied the transport process in the two dimensional electron gas (2DEG) in presence of 
a magnetic field and a dissipative environment at temperature T. By means of imaginary time series 
functional integral method we calculate the decay rates at finite temperature and in the presence 
of dissipation. We have studied decay rates for wide range of temperatures — - from the thermally 
activated region to very low temperature region where the system decays by quantum tunneling. We 
have shown that dissipation and impurity helps the tunneling. We have also shown that tunneling 
is strongly affected by the magnetic field. We have demonstrated analytical results for all the cases 
mentioned above. 
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I. INTRODUCTION 

The decay of a metastable state received an extensive 
study recently 0, 0] because of its important role played 
in physics, chemistry as well as transport in biomolecules 
The tunneling of two dimensional electron gas 
(2DEG) is a general phenomena and is widely encoun- 
tered in mesoscopic systems 0, 0, S : in semiconductor 
heterostructures 0, ^ , in quantum Hall systems ■ 
Here our focus is on tunneling of 2DEG in the presence 
of magnetic field, impurity potential and in a dissipative 
environment. 

In 2DEG the system can be thought of a collection of non- 
interacting charged particles moving in 2D. So one can 
visualize the system in question as if a charged particle of 
mass M moving in a metastable potential U(q) while cou- 
pled to a heat bath environment which is assumed here as 
a collection of harmonic oscillators. In macroscopic sys- 
tems, the tunneling probability is strongly influenced by 
the interaction with the environment at finite tempera- 
ture [TsLll^ . The problem of tunneling in the presence of 
a coupling to many degrees of freedom such as phonons, 
magnons, gives rise to quantum friction which strongly 
affects the tunnehng rate El Q E 13. But all the 
above studies were in the absence of magnetic field. But 
here we shall study the tunneling process of 2DEG in 
presence of magnetic field. Besides dissipation we have 
special attention to the effects of impurity potential on 
the tunneling rate. 

Here we are actually discussing about the decay of a 
charged particle from a metastable potential well. Now 
at high temperatures the decay is thermally activated 
and the rate follows the famous Arrhenius law: 



fciexp{ 



(1) 



where Vb is the barrier height and fd is the attempt fre- 
quency. As the temperature is lowered the classical es- 
cape rate decreases exponentially and it can decay only 
via quantum tunneling. 

Our method of analysis is based on pure thcrmo- 
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FIG. 1: A metastable potential well. 



dynamic equilibrium method which was pioneered by 
Langcr [l^ . In this approach the quantity of interest 
is the free energy of the metastable system. Because of 
states of lower energy on the other side of the barrier, 
the partition function can only b defined by means of 
an analytic continuation from a stable potential to the 
metastable situation depicted in Fig. 1. The procedure 
of analytic continuation leads to a unique Imaginary part 
of the free energy of the metastable state. This quantity 
is then related to the decay probability of the system. 
The thermal quantum rate of decay is related to ImF by 
the following relation: 



K 



-ImF 
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The method we employ to calculate the decay rate is 
a functional-integral approach which is convenient be- 
cause it can easily take into account the dissipation as a 
nonlocal term. Beside this, we are dealing with a prob- 
lem in a higher dimensional space: the particle moves 
in 2-D, and the effective dimension of the dissipative en- 
vironment is infinite. The most convenient method for 
such problem is the path integral method which we use 
here. The Langer's method [ig for the calculation of the 
quantum decay rate was developed through the work of 
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Miller, Stone, Coleman Callan They use the term 

"bounce Euclidean action" Sb and defined the quantum 
decay rate as 



(3) 



where quantum mechanical prcexponential factor ujq is 
related with fluctuations about the bounce. Dissipation 
was first incorporated into the bounce technique by 
Caldeira and Leggett Also Grabert et al |T Hll[l5| 
and Larkin et al [Toj used this method as an effective 
scheme to calculate decay rates in the entire range of 
temperature. In the present paper we extend these 
study by adding a couple of new features in addition to 
the usual dissipative effect: (a) the effect of the magnetic 
field, (b) the effect of the impurity potential. With the 
preceding background I organize the paper as follows. 
In Sec. II I discuss about the free energy and Euclidean 
action of the system. We use the path integral method 
to reduce the infinite dimensional problem to an effective 
one dimensional problem. In Sec. Ill I discuss about 
the thermodynamic method for the calculation of decay 
rates in the high temperature region. In Sec. IV I 
talk about the quantum correction near the crossover 
region. Section V is our main focus. Below the crossover 
temperature the metastable state predominantly decay 
via quantum tunneling. Here we discuss about the effect 
of dissipation, impurity potential and magnetic field on 
the quantum tunneling rate of 2DEG in great details. 
Finally I smimiarizc and conclude the paper in Sec.VI. 

II. FREE ENERGY OF THE DAMPED 2DEG 

I am specifically dealing with two dimensional elec- 
tron gas (2 DEC) which confined in the x-y plane. The 
magnetic field is applied perpendicular to the x-y plane. 
Now following Caldeira and Leggett P, 01 the system is 
assumed to be coupled linearly to its environment which 
can be represented as a collection of harmonic oscillators. 
Following Feynmann- Vernon [20| one can integrates out 
all the environmental degrees of freedom and other ir- 
relevant degrees of freedom. So the problem effectively 
comes out to be 1-D problem. Our starting point is the 
total Hamiltonian of the Global system i.e. the system 
and the reservoir which are coupled linearly: 



H = Hs + Ti-B + Hi, 



(4) 



where subscripts S, B, and / stand for 'system', 'bath' 
and 'interaction' term respectively. 
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where qi = (q^^qy) = q are the co-ordinate of charged 
particle and qj^ = {qj^xiqj,y) are environmental degrees 
of freedom. Now all quantities characterizing the envi- 
ronment may be expressed in terms of the spectral den- 
sity of bath oscillators 



d{uj-UJj). 



(6) 



Sometimes it is needed to introduce cut-off frequency u>c 
0,13. Thus 



Ji{uj) = riiuj-'"exp{ ). 



(7) 



In the present paper we follow the second definition of 
spectral density. The partition function of the charged 
Brownian particle can be written as a functional inte- 
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weighted according to the Euclidean action [l]. lll . lia . l2ll | 



gral |l8l | over periodic paths where path probability is 
eij 

5^[g(r),gir)] = 
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The above action can be transformed into an effective 
action 



S^W{T),q'{r)] = / dr 
Jo 



M ■ 
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dr / dT'/C(r - r')g(r)g(T'), 



where (3 = -jt^- The first term in Eq. (9) represents the 
conservative or reversible motion of the particle whereas 
the second term introduces the dissipation. The damping 
kernel IC{t) may be expressed Fourier series 



/C(r) 
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K.{vn)exp{iUnT), 



(10) 



where Matsubara frequencies 
Fourier coefficients are given by /C(t'n) = ^ 
The partition function is defined as 



Z = I D[q]exp - ^^^fe] 



h 



(11) 
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where the functional integral is over all periodic paths 
with period hp. The free energy is given by 

F = -hnZ. (12) 

P 



III. THERMALLY ACTIVATED DECAY 

Temperature dependence of quantum decay rates in 
dissi pat ive systems was first studied by Grabert et al 
[l^ 0, 0|- They define crossover temperature To at 
which the transition between thermal hopping and quan- 
tum tunneling occurs. For temperature above Tq, the 
main contribution in the hmctional integral (11) arises 
from the vicinity of the time-independent trajectories 
g(r) = 0, where the particle sits on top of the potential 
barrier of the inverted potential, and g(r) = qt,, where 
it sits at the bottom of the well as shown in Fig. 2. a 
periodic path near g(r) = can be expanded in terms of 
Matsubara frequencies {vn) as follows 



+ 00 



x'{t) = ^ X'„exp{ivnT) 



(13) 



Now the action becomes 



S[Y,[] = h(3Uh 



Mh(3 
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where, X''^^ = iy'^-uj^ + 2ujciyn + Wn\'j{\vn\), with frequency 
LUb ~ [tt^]^- Since Aq = hence the integral over 

the amplitude Yq diverges. This leads to an imaginary 
part of free energy. The imaginary part of the free energy 
is given by 



ImF 
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Affleck [221 has shown that above Tq the decay rate given 
by Eq. (2) is modified as follows 



(19) 




Now using the fact that 
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and 
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combining Eq. (17) with Eq. (19) we obtain 
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-uj^ and 
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(20) 

where we made use the definition of crossover tempera- 
ture To = 0, . ujR is the largest positive root of 
the equation wfj -I- ujR{'^{uJii) + 2wc) = w^. Now at high 
temperature T » To the product term in the bracket 
(...) becomes unity, thus we obtain 



FIG. 2: The inverted potential - U(q). 



When we put Eq. (13) into Eq. (9) one obtains 



(14) 



n— — 00 



where, A° = i^f^ + ujq + 2ijJcVn + |i^„|7(|fri|), with weU 
frequency = ^ -^^^ ] ^ and cyclotron frequency lOc — 
1^. Now the partition function Zq can be obtained by 
performing the Gaussian integrals over the amplitudes 
Xn- A periodic path around qij) ~ qb can be written as 



y'iT) 



E Kexp{ii^7iT). 



(15) 
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which is the famous Arrhcnius law. The influence of dis- 
sipation and magnetic field on the classical rate arises 
through memory rcnormalizcd barrier frequency lor. 
As the temperature is lowered, the thermally activated 
decay rate is modified by the quantum correction factor. 
Quantum fiuctuations enhances the decay probability be- 
cause they increase the average energy of the Brownian 
charged particle. Now following Hanggi et al |23| we ob- 
tain the resulting leading quantum corrections are found 
to be given by the formula 



exp 



24 ikeTy 



(22) 



This fq is independent of dissipative and magnetic field 
mechanism. Let us now consider frequency independent 
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damping case i.e. ohmic damping case. For the ohmic 
damping case the product in Eq. (20 ) can be evaluated 
in terms of F functions as follows 



F(1-^)F(1-^) 
F(1-^)F(1-^)^ 



where 



2TTkBT 



and 



(23) 
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For strongly damped system 7 >> woj'^b and T >> Tq 
the ohmic quantum correction factor Eq. (23) becomes 



exp< 



'To 
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(1)]}, (25) 



where \E'(x) is the digamma function and a — On 
the other hand, for intermediate temperatures Tq << 
T « Au^Tq we obtain the quantum correction factor 



(26) 



Now in the third case where 7 << woji^h we have [iJ 
iinhi 



2k bT 



UJo 



exp[Da + Oia^)], (27) 



where a = ^jr^"'' and 
lull, 



vI.(l + ^) + vI/(l_^)_vI.(l+^^)_vl.(l_^^) 

V V V V 
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IV. NEAR CROSSOVER REGION 

In the above section we have seen that the quantum 
corrections become very much important as we approach 
the crossover temperature Tq. In the vicinity of To the 
above mentioned semiclassical approximation fails be- 



cause the eigenvalue Aj = A_ 



vanishes for T ~ Tq. Hence the Gaussian integral over 
the amplitude Fi,y_i diverges. Now to regularize this 
divergent integral I add the higher order terms in ampli- 
tudes Yi, F_i [13 ■ I expand the potential U(q) about the 
barrier top ?7(g) = Ub-lMtuKq-qby + \Mck{q- 
qb)^ ■ Now following Grabert et al I obtain the quantum 
correction factor in the crossover region as follows: 



(29) 



where 



fR 



{2i^c 



Al = 
and B = 



TrM0 \ ■ 
2B I 



f 7(i/„))^'n, 



A more detailed analysis on the crossover region has 
been discussed by many others 0, T?, , although 
they did not consider the magnetic field effect which I 
have considered here. Now I switch over to our main 
focus i.e. the quantum tunneling region. 

V. QUANTUM TUNNELING REGION 

Below the crossover temperature Tq, the tunneling 
through the barrier becomes more probable and the de- 
cay rate is mainly given by the tunneling rate. In this 
region the most probable escape path is the bounce tra- 
jectory qsir). For temperature below To the imaginary 
part of the free energy comes from the so-called bounce 
trajectory [23 ■ The decay rate is given by [l^lT^ITsI l: 



F = 



2Trh 



Do 
D' 



exp 



h J' 



(30) 



where Sb is the action (9) evaluated along the bounce 
trajectory qsir), So is the zero- mode normalization fac- 
tor 

rfi/3 

(31) 



So = M / dTciB^T), 



Do is the product of eigenvalues of fluctuation modes 
about the mctastable minimum, and D'^ is the product 
of eigenvalues of fluctuation modes about the bounce tra- 
jectory with the zero eigenvalue omitted. 
In my case I am interested in the 2D electron gas tunnel- 
ing out of the metastable state g^; = in the x direction. 
The electron gas is moving in the x-y plane and mag- 
netic field is applied perpendicular to the x-y plane. The 
2DEG is coupled with a dissipative environment, so far, 
we have considered arbitrary forms of the metastable po- 
tential. Now I consider the practically important case 

of a cubic potential, Ui{qx) = \Mu}'^(^q^ — which is 

relevant for SQUID or a current-biased Josephson junc- 
tion. I also approximate the impurity potential in the y 
direction by the harmonic potential. Thus U{qx,qy) = 

\Mujl{q^ ~ 90 ) h^'^^yly- convenience I use 

the notation x = {x,y) instead of = {qx,qy) 
tunneling is described by the Euclidean action 



The 



dr 



M{x^ + y^) + iMojcXy + Ui{x) 



MLoy 



^ L2to,- 



(32) 
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where I have used Landau gauge. After the tunneling 
the other degrees of freedoms, y and qj can take any 
allowed values. Hence one can easily integrate out the 
environmental degrees of freedom qj and y coordinates. 
Before doing all these integration over qj and y coordi- 
nates, we perform a Fourier transformation in the time 
interval [0, hP]: 

+00 

(x(r),y(r)) = j-g i^n,yn)e'''-^ (33) 



n— — 00 



1] 



hp 



(34) 



n— — 00 



Now the action can be rewritten in terms of Fourier com- 
ponents as follows: 



Mujl)yny-n 



^ +00 
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nijuj]' 



with Ui{xn) = hp dTU{x{T). After doing integration 
over the environmental degrees of freedom qj n, I obtain 



n— — 00 
+00 



1 ' ^ ]^ 

-j^ 7;{Mvl + Mujl + ^y,n)yny-n 



hp 



+ Ul{Xn) + ^ Y M^cXny-nV 



n— — 00 



with ^i,„ = i/o°°rfw^^^;|^; I = x,y. Here we 
are considering tunneling along x-direction, hence I can 
integrate out the y variables. After doing the integration 
over y variables I obtain the effective one dimensional 
action as follows: 



1 1 

5^ = ^ -{M,^l+^,^n)XnX-n+Ui{Xn){35) 



n— — oo 
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XfiX- 



The above effective ID action can be rewritten as shown 



below: 
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effl 
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,{t)] = / dT 
JO 

dT / dr' 



-Mx\T) + U^{x{r)) 



1 











(36) 

/C(t - r') + .g(r - t') 



(x(r)-x(r')) , 
where the normal damping kernel is given by 



/C(r) 



dujJ{uj) 



cosh[uj{^ - t)] 



sm/i[^^Y^] 



and the anomalous damping kernel is given by 

hp t 



9ir) 



(37) 



(38) 



E 



n—^QO 

,2 



n— — oo 



The first term is tlie delta function and has no contri- 
bution to the semiclassical action, thus finally I obtain 



+ CX3 



(39) 

To derive the normal damping kernel I use the identity 
t^ cosh[uC4-r)] _ Thus we have ob- 



V — ^ 



2 sm/i[^^ J 

tained an effective one dimensional problem. The physics 
due to the influence of the dissipative environment and 
magnetic field is underlined in the normal and anomalous 
damping kernel respectively. Now I shall discuss different 
limiting cases. 

(i)Impurity zero and no dissipation 
In this limiting condition the normal damping kernel 
/C(t) = and the anomalous damping kernel g{T) ~ 
^^^'^"^ Thus the effective Euclidean action becomes 



fi/3 



aE _ 
^eff - 



hp 



dT 



]^Mx^{t) + Ui{x{t)) 



(40) 



hp 



hp 



dT 



hp 



dT'[x{T) ~x{t')] 



'M2 



S^jj depends on magnetic field through ujc and this de- 
pendence is proportional to . Thus for a strong mag- 
netic field the tunneling rate [Eq. (30)] vanishes at very 
low temperature. At zero temperature the effective ac- 
tion Eq. (40) can be written as 



hp 



dT 



\mx\t) + [/i(.t(t)) + \mloIx^{t) 



(41) 



6 



Equation (41) shows that the magnetic field renormahzes 
the potential Ui{x{t) such that the local minima at a; = 
is now more stable. Because the renormalized potential 
is now a double well potential and thus the second local 
minima may be lower than at a; = 0. 

(ii)Finite impurity and no dissipation 

Since there is no dissipation, hence the normal damp- 
ing kernel IC{t) = and the anomalous damping kernel 
is given by 



oE 



dr 



Mx^{T) + Ui{xiT)) 



V 



2tt 

Finally I obtain 
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Thus the effective euclidean action becomes 



aE 



hp 



dr 



-Mx^{t) + Ui{x{t)) 



- — {Mujcf 



hp 



dr 



oE 
^eff 



dr 

Veff 
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where rjeff = t] 



Mx^ 



dr 



+ U,ix{r)) 
1 

dr'-. 



(43) 

[xir)~xir')]\ 



From the Eq. (43) it is well 
understood that the effect of dissipation is to suppress 
the effect of magnetic field. Because of the factor 
we should get finite tunneling at any amount of mag- 
netic field strength. So dissipative environment helps in 
tunneling of 2DEG in presence of magnetic field. 
(i)Finite impurity and finite dissipation 
Since both impurity and dissipation helps in tunneling 
2DEG in presence of magnetic field, so the total effect 
of them will be the same. The effective action is 



ih[ 



UJyhfl^-, 

2 1 



[x(r)-x(r')P 



Now one can easily be able to understand that the tun- 
neling rate is finite for any value of magnetic field. Thus 
the impurity acts in opposite way as that of magnetic 
field. So impurity suppresses the effect of magnetic field 
and it actually enhances the tunneling rate. 

(i)Impurity zero and finite dissipation 

I set the impurity potential i.e. ujy = 0. We know 
that the quantum tunneling predominates at very low 
temperature, i.e. at very large imaginary time limit 
(hp —^ oo). So we have to take large time limit val- 
ues of the normal and anomalous damping kernel 0, |^ . 
In the large r limit the normal damping kernel is given 
by 0, 01 ^(''') = 7?777TT and the anomalous damp- 
ing kernel becomes g(r) = ^^l^M ^ r^-s+i , where 77' = 
2M-f /o°° dz^^^. Thus the effective action for the ohmic 
damping (s=l) case becomes 



■ oE 
^eff 



hp 



dr 



rl 



Mx^{T) + Ui{xiT)) 
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+ 2 
[x{t) 



hp 



dr 

x{r')] 



hp 



rfT'[/C(|T-T'|)+.g(r-T')] 







Now in case of strong impurity the anomalous damp- 
ing kernel is negligible compared to the normal damping 
kernel, because in this case the anomalous superohmic 
kernel decays rapidly compare to ohmic normal damping 
kernel. Hence the effective action becomes: 



oE 



hp 



dr 





V_ 

'27r 



Imx\t) + Ui{x{t)) 



dr 



dr'-. 



r'|2 



(44) 
x{t) -x{t')]'^. 



Equation (44) suggests that dirty sample enhances the 
tunneling rate of 2DEG in presence of magnetic field. 
On the other hand for strong damping case one can ne- 
glect the superohmic normal damping kernel compared 
to ohmic anomalous damping kernel. Thus the effective 
action now reads as: 



oE 



dr 



1 

2 „„ 



Mx^{t) + Ui{x{t)) 



dr' 
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TT 
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|(r-r')| 
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dz 
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Mx^{t) + Ui{x{t)) 



dr 



dr'-. 



r'|2 



(45) 

[x{r)-x{r')r. 



[x(t) — x{t')]'^ From the above expression it is well understood that dis- 
sipation opposes the magnetic field effect. 
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(i)Very strong magnetic field 

Under very strong magnetic field the kinetic energy of 
the electron gas freezes. Thus one obtain the effective 
action as follows 



hp 



dr 



iMuj,{x{T)y[T)) + t/i(x(r)) + ^^Mu^y 



1 



dr 



dT'JC{T~T')[x{T)-x{T')f 



Following Jain ei a/ Q I obtain 
oE _ 



I I dxiMLiJc{i{T)y{T)) 





dr 











r — T 



-2[Xc/(t) -Xci{t')Y 



Finally I obtain 

Scl 



MlOc i dx 
Jo 

+ - dr dr' 

2 ./n ./n IT - r 



2C/i(a;) 



(46) 



7^[a;c;(T) ~Xci{t )Y 



Thus as we increase the impurity potential along 
y-direction Sd decreases and the tunneling rate of the 
2DEG along x-direction increases. 



have provided a complete description of the decay of a 
metastable state extending from the classical regime to 
quantum region. I have covered a wide range of tempera- 
ture region — extending from the high temperature ther- 
mally activated region to a very low temperature quan- 
tum tunneling regime. It has been seen that well above 
the crossover temperature Tq, the decay rate follows the 
famous Arrhenius law. Here the preexponential factor is 
affected by the damping and magnetic field through the 
renormalized frequency tuji. Now as we approach to the 
crossover temperature the quantum correction factor to 
the decay rate becomes very much important. Now well 
below the crossover temperature the system decays pre- 
dominantly by quantum tunneling process. I have shown 
that at zero temperature the tunneling rate of 2DEG 
from the metastable state vanishes in a strong magnetic 
field and in absence of both dissipation and impurity. 
The effect of impurity potential and the dissipation have 
been discussed in great details in Sec. IV. Both the im- 
purity and dissipation subdue the effect of magnetic field 
and enhances the tunneling rate. I have discussed this 
tunneling rate in context of a very relevant cubic poten- 
tial. So my results are also applicable to the physically 
interesting problem of quantum tunneling in SQUID or 
current-biased Josephson junctions. 
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